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Neutral excitons in strained axially symmetric In(Ga)As/GaAs quantum dots with ringlike shape 
are investigated. Similar to experimental self-assembled quantum rings, the analyzed quantum 
dots have volcano-like shapes. The continuum mechanical model is employed to determine the 
strain distribution, and the single-band envelope function approach is adopted to compute the 
electron states. The hole states are determined by the axially symmetric multiband Luttinger-Kohn 
Hamiltonian, and the exciton states are obtained from an exact diagonalization. We found that the 
presence of the inner layer covering the ring opening enhances the excitonic Aharonov-Bohm (AB) 
oscillations. The reason is that the hole becomes mainly localized in the inner part of the quantum 
dot due to strain, whereas the electron resides mainly inside the ring-shaped rim. Interestingly, larger 
AB oscillations are found in the analyzed quantum dot than in a fully opened quantum ring of the 
same width. Comparison with the unstrained ring-like quantum dot shows that the amplitude of the 
excitonic Aharonov-Bohm oscillations are almost doubled in the presence of strain. The computed 
oscillations of the exciton energy levels are comparable in magnitude to the oscillations measured 
in recent experiments. 
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I. INTRODUCTION 

Semiconductor quantum dots (QD) have been widely 
investigated in the last decade because of their potential 
applications for electronic^, photonic^,, spintronic£, and 
quantum computing devices^. Different fabrication pro- 
cedures are adopted to grow these structures, and quan- 
tum dots of various shapes and dimensions have been 
realized. A peculiar example is a quantum ring (QR)££, 
which has a doubly connected topology 7 -. Hence, the 
Aharonov-Bohm (AB) effect^, which arises from a change 
of phase of the particle wave function with magnetic field, 
takes place in the scattering-free limit. In circular ring- 
like structures the AB effect is manifested by the cross- 
ings of energy levels of different angular momenta in the 
ground state, which arc referred to as angular momentum 
transitions, and which lead to oscillations of the electron 
ground energy level with magnetic field 7 -. 

In addition to the AB effect for single particle states, 
the ground energy level of the neutral exciton may exhibit 
oscillations when the external magnetic field through the 
QR varies. This excitonic AB effect was predicted for 
concentric one-dimensional (ID) rings^, where the elec- 
tron and hole are found to accumulate different phases in 
each revolution in the presence of an external magnetic 
field. Also, this simple model predicted that the bright 
exciton state can be turned into a dark one. These bright- 
to-dark exciton transitions produce peculiar oscillations 
of the emission intensity with magnetic field, which are 
referred to as the optical excitonic AB effect^. In a 
three-dimensional (3D) geometry, however, the overlap 
between the electron and hole wave functions arc finite, 



thus the emission intensity does not exhibit periodic os- 
cillations as in ID quantum rings. Rather, due to a small 
polarization of the 3D exciton, small oscillations of the 
exciton ground energy level, which are of the order of a 
fraction of meV, were experimentally observe d 10 ' 11 . 

Most works on the excitonic AB effect has been fo- 
cused on In(Ga)As/GaAs quantum rings, which are self- 
assembled by means of the Stranski-Krastanov (SK) 
growth^ and are strained because of lattice mismatch 
between In(Ga)As and GaAs. Strain is involved in both 
the initial phase of the self-assembly of lens-shaped quan- 
tum dots, and the subsequent dewetting procedure which 
leads to an outward diffusion of In adatoms from the cen- 
tral part of the InAs islands. As a result, the formed 
quantum rings have volcano-like shapes, and exhibit 
substantial compositional intermixing between the dot 
and the matrix^. The redistribution process has been 
demonstrated to be strongly temperature dependent^, 
and QR's of various composition and morphology have 
been fabricated^ - - — . For example, self-assembled quan- 
tum rings exhibiting considerable in-plane anisotropy 
have been explored in Refs. [T3 - [l6l More recently, al- 
most axially symmetric In(Ga)As/GaAs quantum rings 
were fabricated and analyzed^, but no height versus ra- 
dius dependence of these rings was determined. 

When the QR opening is covered by the material of 
the ring, the AB oscillations of the single-particle states 
are considerably modifie d 12 ' 15 ' 17 " —. In view of the re- 
cent experimental discovery of the excitonic AB effect 
in type-I quantum ring o 10 ' 11 , a few issues deserve spe- 
cial attention. First, increasing the height of the QD rim 
with respect to the height of the inner layer [sec Fig. 1] 
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enables larger space for both the electron and hole in 
the rim, and therefore the exciton states could exhibit 
a QR-likc behavio r 14 ' 17 " — . Second, band mixing could 
have an important impact on the hole states, and in 
turn their dependence of the exciton ground state energy 
on the magnetic fiel d 20 ' 22 . Third, recent theoretical re- 
sults have indicated that the spatial variation of strain in 
fully opened quantum rings is beneficial for the appear- 
ance and the magnitude of excitonic AB oscillation o 17 ' 20 . 
Strain could lead to an effective separation of the elec- 
tron and hole, thereby increasing the exciton polariza- 
tion. A similar separation of the electron and hole in 
type-II quantum dots was previously found to promote 
excitonic AB oscillations 2 ^. 

In this paper, we explore how the layer covering the 
quantum ring opening affects the excitonic Aharonov- 
Bohm effect. We assume that the analyzed ringlikc quan- 
tum dot (RLQD) is axially symmetric, as depicted in 
Fig. 1(a). The magnetic field dependence of the electron, 
hole, and exciton states are computed for a few values 
of the rim height and fixed height of the inner layer, as 
illustrated in Fig. 1(b). The strain distribution is ob- 
tained within the approximation of isotropic elasticity in 
the continuum mechanical mode l 24 ' 25 . The electron enve- 
lope functions are computed within the single-band effec- 
tive mass approach, whereas the hole envelope functions 
are extracted from both the axial and spherical approx- 
imations of the multiband Luttinger-Kohn model. The 
exciton states are calculated using an exact numerical 
diagonalization. We will investigate how the size of the 
rim affects the excitonic AB effect. Furthermore, we ex- 
plore effects of band mixing on the magnetic field depen- 
dence of the exciton energy levels. We especially analyze 
how strain affects the mixing of the heavy-hole (HH) and 
light-hole (LH) states and what are the consequences of 
this mixing on the exciton states. 

The paper is organized as follows. In Sec. II, we 
present the theoretical framework to compute the strain 
distribution, the single-particle, and the exciton states. 
The results of our numerical calculations are presented 
and discussed in Sec. III. Our conclusions are given in 
Sec. IV. 



II. THEORETICAL MODELS 
A. Model of mechanical strain 

The analyzed RLQD is composed of In(Ga)As, which 
is embedded in an infinite GaAs matrix, and is strained 
due to the lattice mismatch between In(Ga)As and GaAs. 
An important theoretical model that allows one to cal- 
culate strain was introduced by Dawnes 2 ^, who based 
his derivations on the Eshelby inclusion theor y2£. The 
model was subsequently applied to quantum dots of dif- 
ferent shape and compositio n 24 ' 27 . In this approach, the 
clastic behavior of the material is assumed to be linear 
and isotropic. Also, the dot is assumed to be embedded 
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FIG. 1: (Color online.) (a) Profile of the analyzed ring- 
like quantum dot (RLQD). (b) The radial dependence of the 
RLQD height h for three values of the parameter n that ap- 
pears in Eq. (5) . The cross section of the quantum ring defined 
by Eq. (6) is shown by the dash-dotted line. 



in an infinite matrix, with a lattice constant mismatch 

[a d - a m ) 
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between the dot and the matrix. Here, a m denotes the 
lattice constant of the matrix material, and ad is the 
lattice constant in the dot. The following equation is 
obtained for axially symmetric quantum dots 
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Here, £d(r) is equal to e in the dot and is zero in the 
matrix, \ denotes the Lame displacement potential and 
v is the Poisson ratio. The displacement vector u = 
gradx = §^ e p + §J e z an d the components of the strain 
tensor 
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are determined from the computed x = x(/°, 0! z )i which 
is axially symmetric in the analyzed RLQD. To obtain 
the total solution for the components of the strain tensor 
from Eqs. (2) and (3), the term — £d(r), that accounts 
for the initial compression, must be added to the tensile 
strains, while no correction is needed for the shear strain 
components. Furthermore, the hydrostatic strain in an 
arbitrary shaped quantum dot is given by2^ 

l-2v 

£h y d = £ PP + £00 + £zz = —2 — — £d(r), (4) 

and therefore, for a homogenous material it is constant 
inside the dot and equals zero in the matrix. To avoid 
the division by p in Eq. (3b), is calculated by using 
Eq. (4). 



B. Electron states in the conduction band 

Our aim is to analyze nearly cylindrically symmetric 
quantum rings which were experimentally explored in 
Ref. [ill But their dimensions were not accurately de- 
termined, and therefore we used the cross sections of the 
quantum rings explored in Refs. [l2| and [3 m the planes 
(HO), (100), and (110) [see Fig. 1(b)], and revolved them 
around the [001] direction to produce axially symmetric 
shapes. The volcano-like shape of such a formed RLQD 
is shown in Fig. 1(a), and following Ref. [i~2l it is described 
by the function 



Hp) 
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p < R 
p> R 
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Here, R is the radial position of the rim top, ho is the 
height of the RLQD at p =0, h^ is the wetting layer 
thickness, hu is the reference value for the height of the 
RLQD rim, whereas the parameters 70 and deter- 
mine the inner and outer slope of the rim, respectively 
The dimensionlcss parameter 77 is the ratio of the RLQD 
height and hu- For 77=0.8, 77=1.0, and 77=1.2 the vari- 
ation of h with 77 corresponds to the (110), (100), and 
(110) cross section of the ring of Ref. [HI, respectively. 
The values of the parameters used to generate the curves 
in Fig. 1(b) are: R = 11.5 nm, ho~ 1.6 nm, hoc = 0.4 
nm, hM = 3.6 nm, 70 = 3 nm, and 700= 5 nm. Note 
that the analyzed RLQD is assumed to be axially sym- 
metric, which is justified by the fact that the very exis- 
tence of the AB effect is related to the topology of the 
structure. The deviations due to the in-plane anisotropy 
of the structure shape may only affect the period of the 
AB oscillations^. Furthermore, the difference of the rim 
height for three cross sections shown in Fig. 1(b) is small, 



and therefore one might infer that the AB oscillations of 
the single-particle energy spectra for the three values of 
?7 exhibit small mutual differences. Furthermore, the ad- 
vantage of adopting the axially symmetric model is its 
conceptual simplicity. 

We will compare the electronic structure of an RLQD 
with the one of a fully opened ring, whose cross section 
is described by 




p-R 



p < R- a 
R- a < p < R 
R + a < p 



(6) 

where a is proportional to the full-width at half- 
maximum [see dash-dotted curve in Fig. 1(b)]. 

The electron states are extracted from the single-band 
effective-mass Hamiltonian 



He = (P + qA e ) 



1 



2m* (r e ) 



(p + qA e )+H Ze + V e (r e ). (7) 



Here, m* is the electron effective mass, q is the elemen- 
tary charge, p is the canonical momentum operator, and 
A e is the magnetic vector potential which in the Coulomb 
symmetric gauge has the form A e = B x r e /2. Here, B 
denotes the external magnetic field, which is assumed to 
be uniform and directed along the z axis. The Zeeman 
term Hz e has the form 



(8) 



where ps = qh/(2mo) is the Bohr magneton, mo is the 
free-electron mass, g* is the effective Lande g-factor, and 
<t is the Pauli spin matrix. The effective confinement 
potential V e is given by2£ 



V e (v) = a c {s pp + + e zz ) + V ff, e {r) 
= a c e hyd (r) + V off , e (r), 



(9) 



where a c and V Q // ie are the deformation potential and 
the confining potential due to the offset in the conduc- 
tion band, respectively. Because the hydrostatic strain 
is piecewise constant in the employed model of isotropic 
elasticity, the effective potential of the electron in the 
conduction band shifts rigidly inside the dot, and is zero 
in the matrix. 

Because of axial symmetry, the orbital quantum num- 
ber l e , which represents quantization of the z projection 
of the electron orbital momentum L, is a good quantum 
number for the electron state. Thus the electron envelope 
function is given by 



1 



(10) 



where n e denotes the principal quantum number of the 
electron state. 
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C. Hole states 

The hole states are described by the Luttinger-Kohn 
(LK) Hamiltonian)^ 8 .^ 

HlK = Hlk.Jc + Hlk.e + Voff,h ■ 2^4x4 + &JzHz, (H) 

where HLx,k is the kinetic part of the LK Hamiltonian, 
Hlk,s is the part which describes the influence of me- 
chanical strain on the valence-band electronic structure, 
14x4 is the 4x4 identity matrix, V ff,h is the confining 
potential due to the offset in the valence band, the diago- 
nal matrix J z = diag(+3/2, +1/2, — 1/2, — 3/2) contains 
projections of the angular momentum J (J = 3/2) of the 
zone-center Bloch electrons in the valence band onto the z 
-axis, ft is the Luttinger parameter and W z = Tilo c , where 
lu c = qB I mo denotes the cyclotron frequency. Note that 
effects due to the piezoelectric field are negligible, and are 
therefore discarded in our calculation s 10 ' 12 . The energies 
are measured from the band extrema in the matrix far 
away from the dot boundary, and the energy axis for the 
hole states is directed from the top of the valence band 
downward to the top of the split-off band. 

The kinetic part of the LK model is given b y 28 ' 29 



(12) 
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Here, 71, 72, and 73 denote the position-dependent Lut- 
tinger parameters, 7 = (72 +7 3 )/2, and p, = (73 - 72)/2. 
For B oriented along the z-axis, the magnetic vector po- 
tential for the hole states in the cylindrical coordinates 
is expressed by Ah = ■ Bp^/2. Thus, it is straightfor- 
ward to derive the following expressions for k± and k z : 
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: 1 d + qB Ph \ 

' Phd<t>h 2h J 
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The strain-dependent part of the multiband LK Hamil- 
tonian Hlk.c is given by 
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Similar to BIlk.ei the matrix elements of H s are written 
in cylindrical coordinates: 
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Here, a v , b and d are the deformation potentials in the 
valence band and denotes the biaxial strain. 

We briefly discuss the symmetry of the different parts 
of the LK Hamiltonian in Eq. (jlip . Even though V ff,h 
has axial symmetry, H^x.k and Hlk,c lack axial symme- 
try due to the anisotropy of the bulk electronic structure 
of the constituent materials and the strain anisotropy, 
respectively. The difference between 72 and 73 is usually 
small, and therefore p, is a small number. If p is approx- 
imately taken as zero in Eq. (|13p . the HLK,k becomes 
axially symmetric, which is the well-known axial approx- 
imation of a multiband k ■ p models. In such a way, the 
term responsible for the in-plane anisotropy in the R ma- 
trix element of H^K.k which is proportional to e +2l ^ h is 
removed from the multiband Hamiltonian. The term pro- 
portional to e +2l ^ h also exists in the R e matrix element 
of Hlk,e, an d similar to H^K.k, is responsible for the lack 
of axial symmetry in Hlk.s- Nevertheless, the difference 
\V3b — d\ is usually much smaller than |-\/3& + d|- For ex- 
ample, in Ino.5Gao.5As we have \fib — d = 0.9, whereas 
Vo& + d = — 7.5? 5 . If we approximately take d = v36, 
the term proportional to e +2l ^ h is eliminated, and there- 
fore Hlk,e becomes axially symmetric. Along with the 
axial approximation of H^K,k and the axial symmetry of 
the confining potential, the approximation \/3b — d rj 
makes the LK model fully axially symmetric. 

Therefore, for the axially symmetric RLQD, the appli- 
cation of this model is justified when 172—73 1 -C (72+73), 
\V3b— d\ -C |\/3£>+d|. Such a model does not discard mix- 
ing between different hole states due to the off-diagonal 
matrix elements of H^k^i an d therefore could poten- 
tially be more accurate than our previous approach 2 ^, 
which discarded this influence. The application of the 
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later model was based on the fact that in an axially sym- 
metric quantum dot the off-diagonal terms of the strain 
tensor arc generally small, except close to the dot's lat- 
eral boundar y 25 ! 30 . Also, e pp — w 0, except around 
the top of the rim of the analyzed axially symmetric 
RLQD, therefore the off-diagonal terms of the strain- 
dependent Hamiltonian have a small effect on the effec- 
tive potentials 3 ^. This model has the form 

H L K,eff = H LK> k + V h . eff + kJ z H z , (17) 
which is simpler than Eq. (fTTj) . Here, 

Vh,eff = dia,g(Vhh,eff,Vlh,eff, Vih,eff,Vhh,eff), (18) 
where 

V h h,eff(r) = -a v s h y d (r)--s h y d (r)+^bs zz (r)+V off , h (r), 

(19a) 

b 3 

%, e //(r) = ~a v e hyd (r) + -e hyd {r)--be zz (r)+V off j l (r). 

(19b) 

The splitting between the HH and LH bands is expressed 
by the linear combination of the hydrostatic strain £h yd 
and tensile strain e zz . These expressions for the effective 
potentials of the heavy- and light-hole states will be sub- 
sequently employed to explain how the localization of the 
hole depends on the RLQD geometry and the associated 
preferential direction of lattice relaxation. 

Hlk acts on the multiband envelope function spinor 

<./.*(**) 

r i T 

iT/h xjjh ,j,h ,j,h 

- W n H J zh ,+Z/2i w n h J zh >+l/2> w n h ,f zh -l/2> w n h ,fz h -3/2\ ) 

(20) 

where f z h=lh+3zh is the quantum number of the projec- 
tion of the total angular momentum F=L+J onto the z 
axis, and rih is the principal quantum number. Because 
of the axial symmetry, each envelope function has the 
form 

*$K,t. k jM = ^ llh * h ^ h ,i h , 3 AphiZ h ), (21) 

In addition to the multiband LK model given by 
Eq. (11), we employed the single-band approach which 
relies on the spherical approximation of the multiband 
model. This Hamiltonian has the same form as the 
Hamiltonian for the electron states given by Eq. (7), with 
the subscript e replaced by h, and the sign of charge al- 
tered (q — > —q). The states of the heavy and light holes 
are separately computed in the spherical approximation, 
therefore the respective states of the HH and LH excitons 
are separately determined. When band mixing is taken 
into account, the HH and LH excitons mix and produce 
the multiband exciton. 



D. Exciton states: Multiband approach 

The exciton Hamiltonian has the form 

H x = H e + H h + V c . (22) 

Here, H e and Hh are the single-particle electron and 
hole Hamiltonians, given by Eqs. (7) and (11), respec- 
tively. Vc = —q 2 / (47reeo|r e — r^|) denotes the Coulomb 
potential energy of the interacting electron-hole pair, and 
the relative permittivity e is assumed to correspond to 
In(Ga)As, where the electron and hole are mostly local- 
ized. Furthermore, because the Coulomb interaction is 
spin independent, the electron spin s z of the exciton is 
a good quantum number. As mentioned above the hole 
states are classified according to the z projection of the 
hole total angular momentum f z h, but the Coulomb in- 
teraction mixes different f z h hole states in the exciton. 
However, the motion of the exciton center of mass is ax- 
ially symmetric, therefore the z projection of the exci- 
ton total angular momentum f zx = s z + l e — f z h is the 
other good quantum number of the multiband exciton 
state. For the given s z and f zx , the multiband envelope- 
function spinor of the exciton state is expanded in prod- 
ucts of the electron and hole envelope functions spinors 

= E Z)cj..»..»i..»X.,i.(re)ff^ 1 ., + i a -/..(rfc). 

(23) 

The exciton energy levels are extracted from the secular 
equation 

(Eg + S ne! ; eiSa + E nhiat+ i t -f tx - E x ) 

n e ,n h l c 

(24) 

Here, E g denotes the band gap energy of the matrix semi- 
conductor. We note that the secular equations of the 
I/kbj licit) anQl \fzx — l,n x ,l) exciton states are equal, 
while their energies differ by gl/isB. The energy of the 
lowest optically active (bright) exciton state is conve- 
niently denoted by E Xi ±. For zero magnetic held, due 
to electron spin degeneracy and Kramers degeneracy of 
the hole states, the exciton states are arranged in quadru- 
plets. 

E. Exciton states: Single-band approximation 

Because of the larger effective mass of the heavy hole, 
the exciton ground state is the HH exciton. The single- 
band approximation does not take into account the ef- 
fects of spin-orbit coupling on the valence-band states, 
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and therefore the spins of both the electron and the hole 
are good quantum numbers for the HH exciton. Also, the 
z component of the total orbital momentum l x = l e + lh 
is a good quantum number for the exciton. For the given 
l x and either s z or j z , the exciton envelope function is 
expanded in products of the electron and hole envelope 
functions 



EE 



Cn e ,n h ,l e #n e ,l e (r e ) *n h ,L~l e ( T h) 



(25) 

It is straightforward to show that the spin-dependent 
Zeeman terms contribute to the exciton energy by 
E Z1 = 9Zj^bB/2, where g* a = -sign(s, - jz) 
x (9h ~ (~ l)' s *~^'ffe) i s tne effective exciton g factor. 
Hence, El pir * could be subtracted from the exciton en- 
ergy, which simplifies the calculation. 

The expansion of Eq. (25) leads to the secular equation 
for the HH exciton energy 



{Eg + £^n e ,z c + E nhl i x -i e — E x ) 

x '$le,l' e ' ^ t .< ' $n h ,n' k 

+ E E(^^^)^^-^( r '0l^|^ c , e (r e ) (26) 

n e ,m i e 

x < >I .- la (r h )> = 0. 



For zero magnetic field, because of double spin degener- 
acy of the electron and hole, and the orbital degeneracy of 
the single-particle states, the exciton states are arranged 
in octuplets, except the l x = states which are fourfold 
degenerate. 



F. Exciton radius 

We define the average in-plane exciton radius by 

pf.= / ¥*1( Te ,T h )- (pl + p 2 h - 2 PePh coa(4> e -4> h )) 
Jn e Jn h 

x & x (r e ,r h )dr h dr e . 

(27) 



Here, ^ x (r e , r^) denotes either the exciton envelope func- 
tion spinor in the multiband model or the exciton enve- 
lope function in the single-band approximation. For the 
case of the multiband exciton, inserting the expansion of 



Eq. (23) into Eq. (27) leads to 

P\\ = [ C le,ne,n h ,s z Cl s , n ' c ,n' h ,s z 
Z c n e ,n' e n h ,n' h 

x (l e ,n' e \pl\ l e ,n e ) ■ ^n h ,n' h 

,n e7 nh,s z Cl e7 7i' e ,n' h ,s s 

x (s z + l e - f zx , n' h \p 2 h \ s z + l e - f zx , n h ) ■ S ne ,n' e (28) 

Q e ,n e .rih ,s z Cl e — l,n' e ,n' h , s z (^e — lj I P e I ^ e ' ^e) 

x (s z + l e - f zx - 1, n' h \p h \ s z + l e - f zx ,n h ) 

~~ c te,n e ,n h ,s z c l e +l,n' e ,n' h ,Sz (le + lj n ' e \ P e\ le, n e) 

x (s z +l e - fzx + t,n' h \p h \ s z + l e - fzx,n h )\. 

The expression for py of the HH exciton state is similarly 
derived, 



P\ =E E E l C ^ e ,n h Cl e , K , n ' h 
l c n e ,n' e rih,ri' h 

x (le,n' e \pl\ l e ,n e ) ■ 6~ nh ,n' h 

X (l x / e ,7J^ \ ph \ l x le, n h') ' ^n e ,n' e 
X (Z x - l e + 1, "4 | ~ Z e , "ft) 

(Z e + l,rig |p e | l e ,n e ) 



(29) 



X (Z x - le - 1, "4 |Pfe| - 'e, nfc)]. 

We previously found that pii exhibits oscillations with 
magnetic field whose period is similar to oscillations in 
the exciton energ y 17 i 20 . 

We also define the average single-particle radius and 
vertical position for the given exciton state, 

(Pe(h))= / & x \r e ,r h ) -p e{h) -^ x (r E ,r h )dr h dr e . 



(30) 



and 

{ze{h))=l / & x Hr e ,r h ) ■ z e(h) -^ x (r e ,r h )dr h dr e , 

(31) 

respectively. For the concentric ID ring with the radii 
< p e > and < ph >, the exciton polarization depends 
on < p e > — < P h >£. Furthermore, these two radii 
determine the effective exciton orbital radius \ x , which 
is for the HH exciton defined byi^ 



A, 



m h 



m e / < p h > 2 +m h / < p e > 2 



(32) 



Together with < p e > — < ph > it determines the effec- 
tive surface area (as the one between concentric ID rings) 
threaded by the magnetic flux. Because the hole effective 
mass is larger, mostly reflects the behavior of < p e > 
with magnetic field, as demonstrated in Refs. [13 andflU 
for the case where compositional intermixing between the 
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ring and the matrix is present. Furthermore, by using a 
simple analytical argument^, it could be shown that the 
diamagnetic shift coefficient is dominated by a term 
proportional to X x /fj, e h, where [i e h = m e mh/(m e + m/,). 
The dependence of X x on B describes part of the varia- 
tion of the exciton energy around its average parabolic 
dependence on magnetic field. 



G. Estimation of the oscillations of the exciton 
energy 

To resolve the AB oscillations in the exciton energy, the 
authors of Ref. |H proposed to plot the second derivative 
d 2 E X: i/dB 2 as a function of B. We propose here an al- 
ternative way to estimate the period and magnitude of 
the exciton energy oscillation. For that purpose, we fit 
the energy of the lowest optically active state E Xi % as a 
function of B by a polynomial of the fourth order, 



(E x ,i 



i=0 



iB l 



(33) 



where Cj's are the best- fit parameters. We found that 
the fourth order polynomial describes reasonably well the 
correction from a purely parabolic dependence that is 
found for the ID model. The exciton energy residual, 



SE X , 



Ex 



(E T 



(34) 



is expected to exhibit the AB oscillations. 

An additional figure of merit of the exciton state is the 
in-plane Coulomb potential energy, which is defined by 



V C \\ = — 



47ree p|| 



(35) 



It will be demonstrated that for the lowest-energy bright 
exciton state Veil exhibits oscillations around a linear 
function 



(Vb||> =d + d 1 B, 



(36) 



where dj's are the best- fit parameters. The oscillatory 
behavior of Vc\\ is then resolved by computing the resid- 
ual 



sv. 



c\\ 



v. 



(Vc\\) 



(37) 



H. 



Oscillator strength and photoluminescence 
intensity 



For the fully opened quantum rings, we previously 
found that oscillations of the exciton energy levels oc- 
cur together with oscillations in the oscillator strength 
for exciton recombination 2 ^. For the z-th exciton state, 
this oscillator strength is given by 



fx 



m oE Xi i 



{Uco e 



pKo)| 2 |m| 2 



(38) 



Here, e denotes the unit vector of light polarization, u c q 
and u v q are the periodic parts of the Bloch functions 
of the electron in the conduction and valence band, re- 
spectively, and M denotes the transition matrix element 
between the envelope functions^ 



M = 



n e Jn 



S(r e -r h ) ■ W x (v e ,v h )dr e dv h 



(39) 



We assume that e is oriented in the xy plane (in-plane po- 
larized light), thus the matrix element squared between 
the zone center Bloch states is given by 



{U c0 \£ ■ p| U v0 ) 



<P 2 s 

2h 2 



(40) 



where P denotes the Kane interband matrix element. 

The exciton recombination is strongly polarization 
sensitive. By inspecting the spin part of the Bloch 
functions for the in-plane light polarization, one may 
find that the valence band functions ^ r n , h j. x ,+3/2 and 
&nh,f*x,—i/2 are the spin- up states and therefore can 
optically couple with conduction-band spin-up states, 
whereas &n h ,f M ,+i/2 and &n h ,f M ,-s/2 can optically cou- 
ple with spin-down states in the conduction band. There- 
fore, the spin selection rule for the exciton recombination 
is \s z — j z \ = 1. Furthermore, only equal orbital momenta 
of the electron and the hole give a non-zero contribution 
to the transition matrix element in Eq. (39). Therefore, 
\f zx \ = 1 are the only bright exciton states for the in- 
plane light polarization. Similarly, in the single-band 
model only transitions between states of the same spin 
and ^=0 are optically active. 

Even though our main interest are oscillations of the 
exciton ground-state energy with magnetic field, at finite 
temperature T there exists a finite probability of popula- 
tion of higher exciton states, which could smear out the 
oscillations observed in the ground state. A measurable 
quantity at finite temperature is the photoluminescence 
intensity^!, which we define to be dimensionless, and rep- 
resent it by the oscillator strength thermally averaged 
over all exciton states, 



1 PL 



fx,i exp{-E Xil /k B T) 
£\ exp(-E Xti /k B T) 



(41) 



III. NUMERICAL RESULTS AND DISCUSSION 

We assumed that the RLQD, whose shape is shown in 
Fig. 1(a), is made of Ino.5Gao.5As which is surrounded 
by the GaAs matrix- 1 ^ 12 ' 14 . The values of the param- 
eters of the function which describe the RLQD shape 
are extracted from the measurements in Ref. [l2l . and 
are given in Sec. II. B. For the strain computation, the 
lateral boundary of the wetting layer is assumed to be 
circular with radius 100 nm. The calculations are per- 
formed for 77=0.8, 1.0, and 1.2, which are the cases de- 
picted in Fig. 1(b). As noted in Sec. I, in addition to 
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these RLQD's, we compute the exciton states in a quan- 
tum ring of cross section shown in Fig. 1(b), which is de- 
scribed by Eq. © . This ring is assumed to have a height 
h = Hm =3.6 nm, whereas the inner and outer radii at 
half the ring height amount to i?i=8 nm and i?2=15 nm, 
respectively, thus a = 4.95 nm. 

The values of the lattice constants, the deformation 
potentials, the Luttinger parameters, the energy gaps in 



InAs and GaAs, and the bowing parameters for In(Ga)As 
are all extracted from Ref. Hf| The HH and LH masses 
determined from these values are in good agreement 
with the results from measurements^. Therefore, we 
employed the following set of interpolating formulas to 
compute the values of the Luttinger parameters in the 
In(Ga)As alloy 



J 



In(Ga)As 

7i 



In(Ga)As 
72 



In(Ga)As 

7 3 



1 / 1-x 

2 (^GaA: 
if 1-X 



2 7 GaAs 7 InAs _ ^ 



InAs 



1 



1 — x x 



.GaAs 



4 V7i - - 
If 1-x 



2 7 GaAs 



7 JnAs 



27: 



InAs 



4 \^ 7 p a As _j_ 2 7 GaAs 7 InAs _|_ 2 7 InAs 



72 - aAs )(l-x) 



•0.145a;(l-x) 



0.0202x(l - x) 



+ 0.145a;(l - x) 



0.0202a;(l - x) 



GaAs 



7 MGa)As + w 

(73 nAs -72 nAs )^-0.481x(l- 2 ;). 



(42a) 

(42b) 
(42c) 



Such determined Luttinger parameters are consistent 
with the results for the Gao.47Ino.53As alloy-21. Further- 
more, we note that | v36 — d\ = 0.9, is much smaller than 
|V3b + d| = 7.5, and I72 — 73 1 <C (72+73), thus the use of 
the approximate forms of HLK,k and Hlk,e is justified 
for computing the hole states by the multiband model. 

In our calculations, the effective Lande g factor was 
taken to be equal g* = —0.44, the Luttinger parame- 
ter k = 1.72, which are the values that correspond to the 
GaAs matrij*^, and were found to lead to a better agree- 
ment with experimental measurements^. The heavy- 
hole states were computed by using the value of the effec- 
tive masses as in Ref. [13, whereas the heavy-hole effective 
g factor g^ = —3.56 was selected such that the effective 
exciton g factor g ex = —4, as derived from the recent pho- 
toluminescence measurements on In(Ga)As/GaAs and 
InAs/InP quantum dot a 40 i 41 . The negative exciton g 
factor agrees with the theory of the multiband exciton 
state s 40 ' 42 , and implies that the lowest-energy hole state 
in the single-band approximation is a spin-up HH state. 

The conduction-band offset in the In(Ga)As/GaAs sys- 
tem is assumed to amount to 82% of the band-gap dif- 
ference between In(Ga)As and GaAs*^ The value of the 
relative permittivity is determined from e(x) = 15.1 — 
2.87x + 0.67x 2 (Rcf.[i3), for the assumed molar fraction 
value :e=0.5. The exciton states are computed in the 
range of B from to 30 T. The ambient temperature T= 
4.2 K is assumed. 

In addition to the effective potentials, the electron 
and hole localization could be related to the band edges, 



which are defined by 2 ^ 

E cb {r) = V e (r) = a c e hyd (r) + V f f , e (r), (43a) 

E hh (r) = V offJl (r) +P £ + sign(Q £ ) 

. (43b) 

X VI Qe I 2 + I S e \ 2 + \Re \\ 

Eih{t) = V off , h (r) + P e - sign(Q 6 ) 

/ ( 43c ) 

X VI Qe I 2 + I Se | 2 + I Re | 2 , 

for the conduction-band electron, the heavy hole, and the 
light hole, respectively. 

The band edges in the conduction, HH, and LH bands 
for rj = 1 are shown in the upper panel of Fig. 2 [(Figs. 
2(a-c)]. These diagrams are generated by using the re- 
sults from calculations for the strain distribution and 
Eqs. (43a)-(43c). Hydrostatic strain is constant inside 
the dot and zero in the matrix, therefore the conduction 
band edge, which is shown in Fig. 2(a), is piecewise con- 
stant. Close to the rim top the dot lateral dimension is 
much smaller than its vertical dimension. Thus, the lat- 
tice is mainly relaxed in the lateral plane, whereas closer 
to the bottom of the RLQD there is no preferential di- 
rection for the lattice relaxation. Also, the quantum well 
for the heavy hole is deeper close to the dot center than 
in the rim. As a consequence, the heavy hole is con- 
fined in a shallower effective potential well at the top of 
the rim than at the bottom [see Fig. 2(b)]. Furthermore, 
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FIG. 2: (Color online) (Upper panel) Contour plots of the band edges in the In(Ga)As/GaAs RLQD for rj = 1: (a) E C b, (b) 
Ehh, and (c) Eih- (Middle panel) The probability densities of the ground single-particle states which correspond to the band 
edges in the upper panel for B = 0: (d) the electron, (e) the heavy hole, and (f) the light hole. (Lower panel) The probability 
densities of the ground single-particle states for B = 0, but with strain discarded: (g) the electron, (h) the heavy hole, and (i) 
the light hole. 



due to strain, the barrier for the LH states is erected in- 
side the dot, and the shallow effective potential well for 
the LH states is formed outside the dot, which is shown 
in Fig. 2(c). One might infer that because of confine- 
ment of the HH and LH states in the different regions of 
the strained structure, band mixing is reduced, which is 
advantageous for the application of the single-band ap- 
proximation when computing the exciton states*^ 

Along with the variations in the band edge, the prob- 
ability densities of the ground electron, HH, and LH 
states for the cases 77 = 1 and B = are displayed in 
Fig. 2. The middle panel [Figs. 2(d-f)] displays the elec- 
tron and hole localization in the strained RLQD, whereas 
the lower panel [Figs. 2(g-i)] show the states in the ab- 
sence of strain. The localization of the electron in the 
ground state does not qualitatively depend on the pres- 
ence of strain, which could be inferred from a comparison 



of Figs. 2(d) and 2(g). However, strain leads to a deeper 
effective potential well for the heavy hole in the crater 
than in the rim, therefore the HH ground state is con- 
fined closer to the RLQD center [compare Figs. 2(d) and 
2(e)], i.e. there is an obvious tendency of the heavy-hole 
probability density in the strained RLQD to leak towards 
the center. On the other hand, the heavy-hole in the un- 
strained structure is mainly localized in the rim, like the 
electron [compare Figs. 2 (g) and 2(h)]. The peculiar 
effective potential that confines the light hole shown in 
Fig. 2(c) favors the LH localization above and below the 
dot, as illustrated in Fig. 2(f). However, when strain is 
not taken into account, the LH state becomes mainly 
localized in the dot, as shown in Fig. 2(i). Different lo- 
calization of the LH states implies that band mixing has 
considerably different effects on the hole states in strained 
and unstrained RLQD's. 
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To demonstrate how the effective potential for the 
heavy hole varies with the height of the RLQD rim, we 
plot in Fig. 3(a) Ehh as a function of p at height z=0.8 
nm, which is half the height of the RLQD crater. It is 
obvious that when the rim height decreases, the potential 
well in the inner layer becomes shallower, and the barrier 
in the rim becomes lower. The probability densities of 
the ground f z = +3/2 states computed by the LK model 
in both the strained and unstrained RLQD's, which are 
displayed in Figs. 3(b) and 3(c), are similar to the prob- 
ability densities of the respective HH states, shown in 
Figs. 2(e) and 2(h), respectively. However, due to the 
peculiar in-plane variation of Ehh shown in Fig. 3(a), the 
probability density of the multiband hole state in the 
strained structure leaks more effectively to the RLQD 
center. Therefore, when the rim height increases the 
hole localization inside the rim is reduced when strain 
is present. On the other hand, because the electron ef- 
fective potential is stepwise, the influence of increasing 
the rim height on the electron localization inside the rim 
is more effective. Therefore, such a difference is advanta- 
geous for increasing the exciton polarization, and in turn 
it could increase the excitonic AB oscillations. 

The single-particle energy levels as function of the 
magnetic field for r\ = 0.8 (dotted red lines), rj = 1.0 
(dashed green lines), and r\ = 1.2 (solid blue lines), are 
shown in Fig. 4. For r\ = 1.2 the rim of the RLQD is high 
enough to establish that the electron is localized mostly 
in the rim. For this case, the orbital momentum transi- 
tions arc at B^.U T, B 2 =17.2 T, and 5 3 =27.3 T. One 
might note that the orbital momentum transitions in ID 
rings occur when an odd multiple of half flux quantum 
threads the ring£. Therefore, they are arranged according 
to Bi/Bi = (2i — 1), i=2,3,. . . , where B\ is the magnetic 
field of the first orbital momentum transition. This re- 
lation is approximate for r\ = 1.2, whereas for smaller 
r) deviation from the ID ring case becomes larger. It is 
ascribed to a larger localization in the inner layer when 
the rim's volume decreases. Consequently, the orbital 
momentum transitions shift towards higher B values. 

Because the effective potential well for the HH states 
is deeper inside the inner layer than in the rim [see 
Fig. 3(a)], the hole is localized closer to the RLQD center 
than the electron. Therefore, the first orbital momentum 
transition for r] = 1.2 is at much larger B than for the 
electron [compare Figs. 4(a) and 4(b)]. For RLQDs with 
a smaller rim, which are the cases r\ = 0.8 and r\ = 1.0, 
the hole is confined in a smaller volume, but the effective 
potential barrier for the heavy hole in the rim lowers due 
to strain. A consequence of the later effect is an increase 
of the average heavy-hole radius. Therefore, the first or- 
bital momentum transition between the heavy-hole states 
shifts towards smaller B values when r\ decreases. 

One might note that the ground HH energy level in 
Fig. 4(b) lowers below the B = value, which is a con- 
sequence of a too large Zeeman splitting in the single- 
band approximation of the hole states. However, mix- 
ing between the hole bands due to the off-diagonal ki- 
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FIG. 3: (Color online) (a) Ehh variation with p for z = 0.8 
nm. (b) The probability density of the hole ground state 
when band mixing is taken into account, (c) The probability 
density of the hole ground state when band mixing is taken 
into account, but strain is discarded. 



netic terms of the LK model opposes effects of the Zee- 
man splitting. Furthermore, the off-diagonal strain- 
dependent terms affects the angular momentum transi- 
tion, as Fig. 4(c) shows. As a matter of fact, the diagrams 
of the electron and hole states in Figs. 4(a) and 4(c) ap- 
pear qualitatively similar. On the other hand, a com- 
parison between Figs. 4(b) and 4(c) demonstrates that 
the single-band and multiband models of the hole states 
exhibit considerable qualitative discrepancy, i.e., band 
mixing has profound effects on the Aharonov-Bohm os- 
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FIG. 4: (Color online) The ground state energy of: (a) the 
electron, (b) the heavy hole within the single-band approxi- 
mation, and (c) the hole with band mixing taken into account 
in the LK model. The dotted red lines denote the case r\ = 0.8, 
the dashed green lines are the states in the RLQD for r\= 1.0, 
and the solid blue lines depict the case r\ — 1.2. The angular 
momentum transitions are labeled by dots, and the orbital 
(angular) quantum numbers of the electron and hole ground 
states are explicitly indicated. 



dilations of the hole ground energy level. Therefore, the 
changes of the hole localization due to the peculiar ef- 
fective potential variation with the rim height shown in 
Fig. 3(a) cannot straightforwardly explain how the single- 
particle energy levels in the LK model are affected by the 
increase of the RLQD height. 

Figure 5 shows how the low exciton energy levels vary 
with magnetic field. The HH exciton energies determined 
by the single-band approximation are shown in the left 
panel, whereas the results of the multiband calculation 
are displayed in the right panel. For both models, the 
presence of small oscillations of the exciton ground-state 
energy around a parabolic function arc found in Fig. 5 
when rj increases. This is similar to experimcnt o 10 ' 11 and 



our previous work on fully opened quantum rings i 17 i 20 
For 77 = 0.8 and r\ = 1.0, the exciton ground energy level 
is optically active (cither l x = or f zx = —1) in the 
whole explored range of B from to 30 T. On the other 
hand, for 77 = 1.2 the l x = and l x = — 1 levels cross each 
other at 19.1 T. This crossing is a consequence of the de- 
creased Coulomb interaction between the electron and 
hole due to increased exciton polarization when the rim 
height increases, which was previously demonstrated in 
Figs. 2 and 3. However, the multiband exciton does not 
exhibit the angular momentum transitions in the ground 
state shown in Fig. 5(f). The difference between the HH 
and multiband exciton ground states shown in Figs. 5(c) 
and 5(f) is related to both band mixing and locations 
of the single-particle orbital (angular) momentum tran- 
sitions. In the single-band approach, the electron orbital 
momentum transitions are misplaced with respect to the 
hole orbital momentum transitions. Therefore, a change 
of the exciton orbital momentum l x = l e + 1% is possible 
when B varies. When band mixing is present, however, 
the orbital/angular momentum transitions in the conduc- 
tion and the valence bands take place at similar magnetic 
field values. Thus, no exciton angular momentum transi- 
tion are found in Fig. 5(f), but the lowest bright and dark 
exciton energy levels approach each other with B, and 
eventually the two states cross each other at B = 39.1 
T. Yet, they almost appear as a doublet, with the energy 
difference not exceeding more than 0.16 meV in the range 
from 30 to 50 T. 

For comparison, l x = and l x = — 1 exciton energy lev- 
els in the fully opened quantum ring of the cross section 
shown in Fig. 1(b) are also displayed in Fig. 5(b). For 
convenience, these levels are shifted down by 45 meV. 
One may notice that the exciton ground state is bright 
(l x = 0) in the whole range of B, which is similar to 
our previous finding for quantum rings of rectangular 
cross section 2 ^. Furthermore, the energy difference be- 
tween the lowest optically active and dark state increases 
with B, thus no crossing between the two is observed 
at even much higher magnetic fields. Also, we show in 
Fig. 5(e) f zx = —1 and f zx = —2 energy levels in the 
unstrained RLQD, which are increased by 125 meV. The 
exciton ground energies of both the fully opened quantum 
ring and the unstrained RLQD vary more slowly than in 
the case of the strained RLQD's. The smaller diamag- 
netic shift is due to smaller mean square of the in-plane 
electron-hole separation^, and therefore it demonstrates 
that the exciton polarization increases due to strain, as 
inferred from Figs. 2 and 3. 

The oscillations of the lowest-energy bright exciton 
level with magnetic field in Fig. 5 are resolved by plot- 
ting d 2 E x _i/dB 2 as a function of Bin Figs. 6(a) and 6(b), 
for the cases of the HH and the multiband exciton, re- 
spectively. Similar variations of the exciton energy resid- 
ual SE x ^i, defined in Eq. (34), with B are displayed in 
Figs. 6(c) and 6(d), for the HH and the multiband exci- 
tons, respectively. One may notice that both d?E Xt i/dB 2 
and 6E Xj i exhibit oscillations of decreasing magnitude 
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FIG. 5: (Color online) Variations of the lowest-energy levels of a few orbital (angular) momenta for the electron spin up with 
magnetic field. (Left panel.) The HH exciton energy levels for: (a) r\ = 0.8 (b) rj — 1.0 (the exciton energy levels decreased by 
45 meV in the fully opened quantum ring are also shown), and (c) rj = 1.2. (Right panel.) The multiband exciton energy levels 
for: (d) r] = 0.8, (e) rj = 1.0 (the exciton energy levels increased by 125 meV in the same RLQD but with strain discarded are 
also shown), and (f) rj = 1.2. 



when rj decreases. This is due to the peculiar strain dis- 
tribution which leads to a barrier lowering inside the rim 
and raising inside the inner layer, as Fig. 3(a) illustrates. 
It in turn leads to a decrease of the exciton polarization. 

The oscillations of the multiband exciton are found to 
be larger. In our previous work we found that oscilla- 
tions of the exciton ground energy level are established 
by means of anticrossings with higher exciton states of 
the same orbital momentum^. These anticrossings ap- 
pear to be more effective when band mixing is taken into 
account. Nevertheless, the amplitudes of both the HH 
and multiband exciton energy level oscillations are of the 
same order of magnitude, which is 0.1 meV. It agrees 
well with recent measurements on self-assembled type-I 
quantum ring a 10 ' 11 , but it is larger than what we previ- 
ously computed for cup-shaped islands^. Furthermore, 
because the magnitude and period of the oscillations for 
different rj arc close to each other, and three values of rj 
correspond to three cross sections of the realistic struc- 
ture, we deduce that the lack of axial symmetry would in- 
deed have a small effect on the excitonic AB variations, as 
inferred in Sec. I. One might deduce that for the case of 
the structure in-plane anisotropy, as one observed in the 
experiment of Ref. UH, the variations of both d 2 E Xy i/dB 2 
and 5E Xi x with B are close to the average height of the 



ring, i.e., the case rj = 1. The quantitative estimation of 
the effects of the in-plane anisotropy of the shape of the 
structure is beyond the scope of the present analysis. 

Also, for the geometry of the quantum ring, shown in 
Fig. 1(b), the oscillations shown by the dash-dotted lines 
in Figs. 6(a) and 6(c), demonstrate that the excitonic AB 
effect in the fully opened case is smaller than the oscil- 
lations in the strained RLQD. Larger oscillations in the 
latter are mainly a result of the increased exciton polar- 
ization as a consequence of the presence of the inner layer 
and the associated peculiar strain distribution. More- 
over, for rj = 1, the unstrained RLQD exhibits smaller 
oscillations than the strained RLQD, as shown by the 
dash-dotted lines in Figs. 6(b) and 6(d). Therefore, com- 
parisons between the unstrained and strained cases in 
Figs. 6(b) and 6(d) indicate that the presence of strain 
is not detrimental for the existence of the excitonic AB 
effect, but it leads to an almost double enlargement of 
this effect. 

To explain in more detail the origin and shape of the 
exciton energy oscillations with magnetic field, the char- 
acteristic parameters of the lowest-energy bright HH ex- 
citon state as function of B are displayed in Fig. 7 for 
77=1.0. We first note that variation of the average elec- 
tron radius < p e > shown in Fig. 7(a) is much smaller 
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FIG. 6: (Color online.) d 2 E x , 1 /dB 2 as function of B for: 
(a) the HH exciton (the dash-dotted magenta line shows the 
result for the fully opened quantum ring), and (b) the multi- 
band exciton (the dash-dotted magenta line shows the result 
for the unstrained RLQD). Variation of the exciton energy 
residual 5E Xt i with B for: (c) the HH exciton and (d) the 
multiband exciton. 




FIG. 7: (Color online) Variations of the parameters of the 
lowest-energy bright HH exciton state with B for r\ — 1.0: (a) 

< p e > (solid red line) and < ph > (dashed blue line) radius, 
(b) < z e > (solid red line) and < Zh > (dashed blue line), (c) 

< p e > — < Ph > (solid magenta line) and < z e > — < Zh > 
(dashed green line), and (d) pu (solid magenta line) and 
(dashed blue line). 



than the variation of the average hole radius < ph >, 
which could be ascribed to the smaller electron effective 
mass. Variations of the vertical electron and hole posi- 
tions with B shown in Fig. 7(b) are practically negligible, 
and also the average vertical positions of the electron 
and hole are almost equal. However, < p e > — < p^ > 
decreases considerably with B, which is demonstrated 
in Fig. 7(c). This decrease is a result of the increas- 
ing Coulomb interaction when magnetic field increases. 
Therefore, the exciton polarization decreases, and conse- 
quently the excitonic AB oscillations abate, as demon- 
strated in Figs. 5 and 6. Similarly to < p e > — < >, 
the average in-plane exciton radius p» exhibits oscilla- 
tions around a decreasing function of B, as shown in 
Fig. 7(d). However, due to the increase of < p e > with B 
and the smaller electron effective mass, the effective exci- 
ton radius X x is an increasing function of B. Furthermore, 
it oscillates similar to py , but its change is much smaller, 
just 0.3 nm when B increases from to 30 T. Also, in this 
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FIG. 8: (Color online) (a) The in-plane Coulomb potential 
Vc\\ (solid blue line) oscillates around (Vc||) which is modeled 
by a linearly decreasing function of B (dashed green line), (b) 
Variation of the residuals of the in-plane Coulomb potential 
8 Vc\\ (solid blue line) and the exciton energy 8E x ,i (dashed 
green line) with magnetic field, (c) The oscillator strength 
for recombination of the lowest-energy bright exciton state 
(solid blue line) exhibits oscillations around (fx,i) the linear 
function of B (dashed green line). The results for the lowest- 
energy bright HH exciton state for 77=1. are displayed. 
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FIG. 9: (Color online) (a) Photoluminescence intensity Ipl 
for recombination of the HH exciton as function of B (the re- 
sult for the fully opened quantum ring is shown by the dash- 
dotted magenta line), (b) Ipl for recombination of the multi- 
band exciton (the result for the unstrained RLQD is shown 
by the dash-dotted magenta line). 



magnetic field range the increase of X x is smaller than the 
decrease of < p e > — < ph >, and thus one may deduce 
that the effective surface area is diminished. 

Because the RLQD's electron and hole states are better 
confined when B increases, the overlap integrals between 
the single-particle states increase with B. It leads to the 
decrease of p\\ shown in Fig. 7(d), and in turn Vc\\ de- 
creases. Also, it exhibits oscillations around the linear 
function of B, as Fig. 8(a) shows. By careful inspection 
of Figs. 7(d) and 8(a), one may find that the oscillations 
of Vc\\ are in phase with the oscillations of pn. The 
excitonic Aharonov-Bohm oscillations are well resolved 
in Fig. 8(b), where the plots of SV C \\(B) and SE xA (B) 
are both shown. These two quantities oscillate oppo- 
site to each other, illustrating the fact that the Coulomb 
interaction tends to suppress the excitonic Aharonov- 
Bohm oscillations. As a consequence, the oscillations 
of 6E Xt i(B) have a smaller amplitude than the oscilla- 
tions of 8Vc\\(B). Figure 8(c) displays that opposite to 
the variation of Vq\\ with B, enlarged confinement of the 
single-particle states causes the oscillator strength for re- 



combination of the lowest-energy bright exciton state f x ^\ 
with B to increase. Furthermore, f x ,i exhibits oscilla- 
tions around a linear function of B which are opposite to 
the oscillations of Vc\\- Those oscillations are nicely cor- 
related with the oscillations of p\\ shown in Fig. 7(d). For 
example, whenever p» exceeds the linear fitting function 
of B, drops below the similar linear dependence on 
B. 

Figures 9(a) and 9(b) show how the photoluminescence 
intensity varies with B for the cases of the HH and multi- 
band exciton, respectively. Because the exciton is dom- 
inantly localized in the bright exciton ground state for 
small B, the Ipl dependence on B is similar to the vari- 
ation of f x ,i shown in Fig. 8(c). However, Ipl decreases 
for 77 = 1.2, because of a finite population of higher ex- 
citon dark states at finite temperature and the smaller 
energy difference with the bright exciton states. Also, all 
cases shown in Fig. 9, including the RLQD where strain 
is discarded, exhibit an oscillatory variation of Ipl with 
B. 
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IV. SUMMARY AND CONCLUSION 

We explored how the geometry affects the neutral exci- 
ton states in quantum dots whose shape resembles rings 
but have a layer inside the nominal ring opening. Our 
calculations reveal that such an inner layer enhances the 
excitonic AB oscillations for the In(Ga)As/GaAs system, 
and are larger than in the case of fully opened quantum 
rings. We show that an increase of the height of the 
quantum-dot rim has the opposite effect on the orbital 
(angular) momentum transitions, which is ascribed to ef- 
fects due to strain. As a matter of fact, we found that 
strain leads to an increased separation between the elec- 
tron and hole, and also it reduces the mixing between 
the HH and LH states. The results of the single band 
approximation are found to compare favorably well with 
the multiband calculations, except that slightly larger os- 
cillations of the exciton energy levels are found if band 



mixing is taken into account. Even a small increase of the 
rim height is found to bring about a considerable shift of 
the angular momentum transition between the exciton 
states. The magnetic field dependence of the photolu- 
minescence intensity is also affected by the presence of 
the inner layer and the variation of the rim height. The 
magnitude of the computed excitonic AB oscillations are 
found to be comparable to those measured experimen- 
tally. 
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